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ABSTRACT 


This report describes the first year's work on NASA Contract 
NASS- 27863 entitled "Study of an Instrument for Sensing Errors in a 
Telescope Wavefront." It contains the partial results of theoretical 
and experimental investigations of different focal plane sensor con- 
figurations for determining the error in a telescope wavefront. 
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CHAPTER I 


INTRODUCTION 

The purpose of our work at the Optical Sciences Center is to 
conduct experimental and analytical research on optical testing concepts 
suitable for application as figure sensors for the Large Space Telescope 
(LST) program. 

An optical testing device for this application, which we will 
refer to for brevity as a figure sensor, should have a wide dynamic range 
from a few waves, so that it may serve as an absolute standard for 
alignment-associated problems, to about a hundredth of a wave, so that it 
may serve as a means of sensing for nulling out figure errors or for 
measuring wavefront errors for post-detection image processing. For 
these purposes either selected wavefront samples or the entire system 
wavefront should be accessible for measurement. In order to accomplish 
these goals the figure sensor has to be placed at the focal plane of the 
telescope, and a real star source has to be used. 

As there is no single test that has the range required for 
alignment-associated wavefront errors and also the sensitivity required 
(1/lOOX) for image processing and figure correction, we have divided our 
work into a study of: (1) a coarse range geometrical sensor operating 

from several waves to a large fraction of a wave, and (2) a fine range 
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sensor useful from a large fraction of a wave to approximately 1/lOOX. 

We have narrowed our study to one candidate in the coarse range (a modi- 
fied Hartmann test) and to three contenders in the fine range region 
(Zemike, shearing polarization test, and Zernike polarization test) . 
These are discussed in detail in this report. We expect that many of the 
components used in the fine range sensor will also be suitable for use 
with the coarse range sensor and that the range of operation for coarse 
and fine sensing will have significant overlap. ' 



CHAPTER II 


ABERRATED TELESCOPE WAVEFRONT SIMULATOR 

1. General. 

In order to experimentally determine the performance capabilities 
of the various types of figure sensory we have designed and fabricated a 
device to produce a simulated aberrated telescope wavefront. This simu- 
lator includes a variable magnitude star, imaged by a specially designed 
telescope which can introduce various alignment, focus and figure errors 
such as may occur in the operational LST. The aberrated star image pro- 
duced can then be evaluated by the figure sensor under study. A sketch 
of the simulator appears in Fig. 1. The star source consists of a 
quartz iodine tungsten filament lamp, which is run by a stabi- 
lized power supply. The filament is focused onto a 25 ym pinhole, 
through a narrow band interference filter if required. The star source 
can be moved in the x,y direction to simulate wavefront tilt in any 
direction. This light is collimated by a 1.2m (48”) focal length tele- 
scope objective to simulate a star at infinity. The next 5 elements (bi- 
convex lens, biconcave lens, two flat plates and thick meniscus) are 
collectively called the aberration generator. This group of elements 
produces an f/6 beam into which various orientations and magni- 
tudes of coma, astigmatism, and spherical aberration may be introduced in 
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Fig. 1. Aberrated telescope wavefront simulator. 
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an exact manner. A relay system then reimages the £/6 image. The first 
element of the relay system is a 40 mm cemented doublet which collimates 
the aberrated star image. This element can be moved along the optical 
axis of the system in order to correct for focus shift produced by the 
aberration generator, or in order to introduce known amounts of defocus 
into the wavefront. The second element of the relay system is an 80 mm 
focal length telescope objective, which converges the light to an f/12 
cone, completing the simulation of an f/12 telescope where tilt, defocus, 
astigmatism, coma, and spherical aberration can be introduced into the 
beam in differential amounts of about 1/lOOX and over the range of ±2 
waves. The final collimating lens is related to the figure sensor under 
test and will be discussed in detail later. 

2. Aberration Generator. 

At the heart of the aberrated telescope wavefront simulator is 
the aberration generator. This instrument has been designed .to introduce 
precisely known amounts and orientations of different primary aberration 
types. What follows is a discussion of the design concept, the fabri- 
cated optics, the optical evaluation and the mechanical assembly of the 
aberration generator. 

The optical design is shown diagrammatical ly in Fig. 2. The 
instrument consists of an achromatic air-spaced doublet with a 25.4 mm 
clear aperture and a 152.4 mm focal length. The air-spaced doublet is 
designed for high coma but minimum spherical and chromatic aberration. 
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The residual astigmatism, inherent in a doublet design, is unavoidable 
but can be readily taken into account. This doublet is tilted about its 
rear nodal point (located 12.90 mm inside the last radius R^) to generate 
approximately 2 waves of third-order coma at a tilt of 1.5°. About 0.4 
waves of sagittal astigmatism are generated unavoidably at this angle. 
However, since the coma varies linearly with tilt, while the astigmatism 
varies quadratically, the astigmatism will be negligible for small 
amounts of coma (i. e., for up to about quarter-wave of coma, the residual 
astigmatism is less than 1/lOOX). Higher order aberration terms are 
negligible for the objective. Although the doublet is color corrected, 
for precise work a 40 to 50 nm bandwidth filter should be used. Fol- 
lowing the doublet are a set of parallel plates which can be counter 
tipped to introduce astigmatism. No third-order coma is introduced, and 
higher order aberrations are negligible. These plates can be tilted up 
to 15° at which point 2.2 waves of sagittal astigmatism are introduced. 
Again, this astigmatism contribution varies quadratically. The last 
element, a thick meniscus lens, can be moved longitudinally for genera- 
tion of +2.8 to -1.2 waves of spherical aberration. However, the effect 
of moving this meniscus is complicated since both fifth- and third-order 
spherical aberration are changed. A change in focal position and effec- 
tive focal length must also be taken into account. 
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5. Design Evaluation, 
a. Coma. 

The air-spaced doublet has been ray traced at a field angle of 
1.5° and the important coefficients for two wavelengths appear below in 
Table I. 


TABLE 1 

DESIGN COEFFICIENTS FOR DOUBLET 


X = 546.1 nm. 

Field Angle 1.5° 





SA3 

SA5 

SA7 

CMA3 

CMA5 

AST3 

-1.84587E-3 

+1.70755E-3 +1. 

86313E-4 

-2.23989E-3 

-1.7334E-5 

-8.95394E-4 

X = 632. 8 nm 






-2.57982E-3 

1.66305E-3 1. 

87093E-4 

-2.22001E-3 

-2.4S111E-5 

-8.92584E-4 


The meniscus element is positioned to minimize any spherical 
aberration residual, since extrinsic higher order aberration contributions 
are small and higher order coma and astigmatism are negligible, we can 
compute the Optical Path Difference (OPD) changes introduced by the 
doublet by using 

°^°Doublet " "T CMA3 i>3cos(|>U - j AST3(cos 24,+ l)r2u^, 
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where r is the fractional pupil height, <|) is the angle measured from a 
direction perpendicular to the tilt axis, and U is the fractional tilt 
angle. Below is the OPD calculated for these tilt effects at two wave- 
lengths. As mentioned before, for OPD < wave, the astigmatism contribu- 
tion is < 1/lOOX and can be ignored. 


b. Astigmatism. 

The two plane-parallel plates were ray traced at an angle of 15° 
in an £/6 convergent beam. Because of the high degree of correction in 
the objective, it is expected that extrinsic contributions due to objec- 
tive imperfection will be negligible. Since two counter- tipped plates 
are used the third-order coma cancels and the third-order astigmatism 
dominates, with higher order aberrations inconsequential. The meniscus 
element is positioned to minimize the residual spherical,. A summary of 
the third- and fifth-order coefficients for the two counter- tipped plates 
appears in Table II. The OPD due to astigmatism is calculated from 

“'’“Astigmatism = -7 


c. Spherical aberration. 

The position of the meniscus element relative to the doublet will 
determine the amount of spherical aberration in the system. This aberra- 
tion contribution is strictly speaking very complex consisting of third, 
fifth, seventh and higher order spherical coefficients to determine the 
OPD in the pupil. However, from comparisons of exact ray trace OPD to 
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TABLE II 

DESIGN COEFFICIENTS FOR ASTIGMATISM 


X = 546.1 nm, Plate 

Tilt 15° 



SA3 

SA5 

AST3 

AST5 

4.68100E-4 

-1.50419E-6 

4.83740E-3 

-1.60639E-6 

X = 632.8 nm 



! 

4.67761E-4 

-1.49766E-6 

4.833908E-3 

-1.599410E-6 


that of OPD calculated by using third-, fifth- and seventh-order coeffi- 
cients, the computational error is always better than 5% or 7% in the 
OPD, where the smaller percentage refers to meniscus position changes 
made to introduce small errors (~1/10X) and the larger percentage refers 
to meniscus position changes made to introduce larger errors (~1/2X). 

The OPD is calculated from 

OPD =--i SASr** - SASr^ - A SAl'r^. 

o IZ lb 

As we reposition the meniscus we also change the back focal position 
(BFP) and the aberration generator effective focal length (GEFL) . This 
change in BFP is compensated for in the relay optics and the change in 
GEFL results in an f# change which can be compensated for by stopping 
down the pupil in order to maintain an f/6 cone in the generator and an 
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f/12 cone in the simulator. Figure 3 illustrates the BFP as a function 
on meniscus position. Figure 4 shows the GEFL as a function of meniscus 
position. Figures 5-18 show the variation of third-, fifth- and seventh- 
order spherical aberration with meniscus position for two wavelengths. 

d. Overall optical path difference . 

The graphical coefficients were obtained by ray tracing at vari- 
ous positions of the meniscus. In order to obtain analytical expressions 
for BFP, GEFL, SA3, SA5, and SA7 each of the curves was fitted to a 
polynomial expression which resulted in the following equations; 

SA3 = (2.4683572 x 10"^)x^ - (5.4403322 >: 10“^)x2 

+ (3.9624001 X 10~^)x - 0.092684628 

SA5 = (5.4406814 x 10"®)x3 - (1.0583634 x 10‘^)x^ 

+ (7.0274608 x 10"‘*)x - 0.014209098 

SA7 = (7.9180215 x 10”^)x^ - (1.4642964 x 10”^)*x^ 

' + (9.3790299 x 10"^)x - (1.9275608 x lO"^) 

BFP = -(3.9558132 x 10"^)x^ - (8.0375794 x 10"‘^)x2 

+ (0.30452283)x + 108.12638 

GEFL = (3.8922327 x I0’^)x^ + (3.9824714 x 10"‘*)x^ 

+ (0.65624124)x + 118.61591 

where x is a term in millimeters, which referes to the spacing between 
the plane-parallel plates and the meniscus element. It is found by 
adding the reading on the dial indicator to a constant C. This constant 






Fig. 4. Aberration generator effective focal length (GEFL) versus meniscus position. 
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Fig. 5. Third-order spherical aberration as a 
function of meniscus position for 
X = 590 nm. 
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Fig. 6. Third-order spherical aberration as a function of meniscus 
position for X = 590 nm. 
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Fig. 7. Third-order spherical aberration as a function of 
meniscus position for A = 590 nm. 





Fig. 8. Fifth-order spherical aberration as a fimction of meniscus position for 
X = 590 nm. . 
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Fig. 9. Fifth-order spherical aberration as 
function of meniscus position for 
X = 590 nm. 
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Fig. 10. Seventh-order spherical aberration as a function of meniscus 
position for X = 590 nra. 




20 





Fig. 11. Seventh-order spherical-aberration as a function of meniscus 
position for X = 590 nm. 



Third order spherical aberration (mm) 






Third order spherical aberration (mm) 




6328 A 


24 



X X 


T- Ln 

(ujuu) uojiBJjaqe leouaqds japjo j 


Fig. 15. Fifth-order spherical aberration as a function of meniscus position for 
X = 652.8 nm. '• 
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is to be determined experimentally in the calibration of the generator, 
so as to compensate for residual manufacturing errors. Using these 
equations and the expressions for OPD contribution of the various com- 
ponents of the aberration generator we can generate an overall expression 
for OPD introduced by the total aberration generator as follows: 


OPD„ ^ , = OPD_. , , + OPD,, + OPD„, 

Total Pinhole Doublet Plates 


+ OPD,, . + OPD„ , 

Meniscus Relay 


where the optical path difference for the total aberration generator 
OPD.j,^^^^ is equal to the sum of the OPD's of the individual components of 
the system. 

The optical path difference contributed by lateral shifts in the 
pinhole position ®P*^pinhole effect of tilting the wavefront. 

This is expressed as : 


OPD 


Pinhole 


R 


6r cos(ij)-a) 


12.7 

R 


3r cos(()>-a) 


where h is the radius of the entrance pupil, R is the focal length of the 
aberration generator, 8 is the lateral distance the pinhole is located 
from the axial position, r is the lateral distance from the optical axis 
as a fraction of the pupil radius for the point on the wavefront being 
evaluated, (p is the angular position of the point on the wavefront as 
measured from the vertical y axis, and a is the angular position of the 
pinhole as measured from the vertical y axis. 3 and a are determined by 
the readings of the micrometers that position the pinhole x-y movement 
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stage. They are: 

CLEFL 

a = arctan ^ = 

where CLEFL is the focal length of the collimating lens, and AX and AY 
are the readings on the micrometers for the X and Y axes respectively, 
and AXO and AYO are the positions of the micrometers when the pinhole is 
located on the optical axis, r and 41 are determined by the readings of 
the micrometers on the x-y stage that positions the sampling aperture. 
The measurements made for r are divided by the radius h' of the pupil in 
the plane of the sampling aperture. 


v/(AX-AXO)^ + (AY- AYO) 2 
CLEFL 


arctan 


AX-AXO 
AY- AYO 


, , (SEFL) (RIEFL) , 

^ ~ (GEFL) (R2EFL) ^ 

where SEFL is the focal length of the sampling lens* RIEFL and R2EFL are 
the focal lengths of the first and second relay lenses respectively, and 
GEFL is the focal length of the aberration generator as determined above. 
Then, 

r = ^ /x2 + y2 = i /(PX-PX0)2 + (PY-PY0)2 


<b = arctan — 

y 


= arctan 


PX-PXO 

PY-PYO 


where PX and PY are the readings of the micrometers that position the 
sampling aperture and PXO and PYO are the readings when the sampling 
aperture is centered about the optical axis. 

*This lens recollimates the light from the second relay lens as in 
Fig. 1. 
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The optical path difference contributed by the tilt of the 
doublet OP^Doublet expressed by 


OPD 


Doublet ' - ’l(™A3)r!'u cos9 


= - %(CMA3)r^(U cosij) + U sincji) 

X y 

= - %(CMA3)r3(-g-^){ (GXO-GX)cos(j) + CGY-GYO)sin<j)} 

where Q1A3 is the amount of coma introduced by tilting the doublet 1.5°, 

U is the angle of tilt from the optical axis as a fraction of the maximum 

angle of 1,5°, and 0 is the angular position of the point being sampled 

on the wavefront as it relates to the direction in which the coma is 

being introduced. U and U • are the x and y-axis components of the angle 

X y 

of the doublet and are measured by GX and GY as the readings on the 
micrometers of the orientation device. GXO and GYO are the micrometer 
readings when there is no tilt in the doublet. The 0.040 factor is the 
micrometer movement for a tilt of 1.5° and so gives the angle as a frac- 
tion of the angle for which CMA3 was determined. 

The optical path difference contributed by the tilting and rota- 
ting of the plane-parallel plates is expressed by 


OPD 


Plates "" ■ ^CAST3)a)2r2cos20» 


J5(AST3) 


2CBZ-BZ0) 
15 ( 


r2cos2[(|) - CRZ-RZO)], 
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where AST3 is the amount of astigmatism introduced when the plates are 
tilted at an angle of 15° from their position of normal to the optical 
axis, 0 ) is the fraction of the 15° that they are actually tilted, and 9' 
is the angular position of the point being sampled on the wavefront as it 
relates to the direction in which the astigmatism is being introduced and 
is equal to the difference between (|) and the angle from the vertical y 
axis the plates are rotated. BZ is the reading on the tilt dial indica- 
tor showing the tilt and BZO is the reading when the plates are both per- 
pendicular to the optical axis. The difference is multiplied by two 
because each unit division on the dial indicator is equal to 2° then 
divided by 15 to show the angle as a fraction of the angle for which AST3 
was determined. RZ is the reading on the rotation indicator, and RZO is 
the reading when the plates tilt so as to introduce sagittal astigmatism 
in the direction of the y axis. 

The optical path difference contributed by the longitudinal move- 
ment of the meniscus lens 0PD„ . is represented by 

Meniscus ^ ^ 

i 6 r2 - (SA3)r‘* - ^ (SA5)r^ - ^ CSA7)r8 

Meniscus 2 Q 2 z 8^"^ 12^'^ 16 

fl2 7")^ 1 

= - 2 ( GEfep - (BFP-BFP0)r2 - 1 (SA3)r- 

- i (SA5)rB - ^ (SA7)r8, 

where 6 is the change in focal position from the paraxial spherical 
aberration condition position with change of the meniscus lens position 
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and Q is the focal length of the aberration generator. BFP is the back 
focal position of the star image and BFPO is the back focal position for 
minimum spherical aberration. SA3, SA5, SA7, GEFL, and BFP are defined 
in detail above. 

The optical path difference contributed by the movement of the 
first relay lens OPD^glay expressed by 

“'’“Relay ° 

fl2 71^ 

2(GEFLy2 (DZ0-DZ)r2, 

where is the movement of the first relay lens. DZ is the reading of 
the dial gauge for the position the relay lens is at and DZO is the lens 
position when its focal point coincides with the focal point of the 
aberration generator with the meniscus lens positioned to give minimum 
spherical aberration. 

4. Mechanical Components. 

The air-spaced doublet in the aberration generator is mounted in 
a cell that threads into a Lansing Research Model 10.203 angular orienta- 
tion device with linear micrometer adjustments. Each 0.0001" division of 
linear motion of the micrometers corresponds to a 13.5 arc second rota- 
tion of the lens. The cell is threaded into the orientation device until 
the rear nodal point of the lens coincides with the intersection of the 
X- and y-gimbal tilt axes. The direction and amount of coma introduced 
is calculated from the angle at which the doublet is tilted from the sys- 
tem optical axis. 
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The two plane-parallel plates are mounted into metal plates on 
individual axes parallel to one another. A toothed gear on each plate 
meshes with the gear on the other plate. Since both gears are of equal 
size, a tilt of one plate will cause an equal and opposite tilt of the 
other plate. One of these plates is geared to a dial indicator. Each 
0.02 division on the dial is equivalent to a 2.4 minute tilt of each 
plate. This whole assembly is bearing-mounted so that it will rotate 
about the system optical axis. The bearing rotation is gear controlled. 
These gears also operate a digit counter where a one digit change corre- 
sponds to 1 degree of rotation about the optical axis. The amount of 
astigmatism is computed from the angle difference between the two plates, 
and the direction by tilt orientation. 

The thick meniscus lens is mounted into a cell that is attached 
to a micrometer slide. Initially the lens is positioned so as to give 
minimum spherical aberration for the system. The lens can then be moved 
toward or away from the air-spaced doublet and parallel plates to put 
various amounts of positive or negative spherical aberration into the 
system. The travel of the lens is measured with a dial gauge having a 
25-mm travel and reading to 0.01 mm. The amount of spherical aberration 
is calculated from the distance between the adjusted position of the 
meniscus lens and its nominal position when the system has minimum 
spherical aberration. 

The two cemented doublets for the relay optics are mounted in 
individual cells. The cell for the first lens is mounted on a micrometer 
slide and is positioned so that its focal point coincides with the image 
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formed by the aberration generator. As the image point moves along the 
optical axis when the meniscus lens is moved, the first relay lens is 
moved on its micrometer slide so as to keep its focal point and the image 
in coincidence. This keeps the beam coming from this lens collimated. 

The collimated beam is then brought to a focus by the second lens, an 
f/12 telescope object. A photograph of the complete assembly appears in 
Fig. 19. 






Fig. 19. Aberration generator and relay system. 


CHAPTER III 


THE COARSE RANGE SENSOR 


1. General. 

Several types of coarse range (geometrical) optical tests were 
considered. Three of the most common are the Foucault knife-edge test, 
the wire test and the Ronchi test. Effectively, in each of these types 
of tests, a particular type of obscuration mask (knife edge, wire, multi- 
ple wire) is placed in the vicinity of the focal plane. Knowing the mask 
parameters and location and by measuring shadow locations on the pupil, 
one can obtain ray slope information and, by integration, wavefront data. 
However, these tests suffer from a common flaw in that the test mask has 
to be rotated for nonsymmetric aberrations like astigmatism and coma. 
Also, to obtain meaningful error measurements we need to scein the pupil. 

A test that is at least as sensitive as the above and which does not 
require mask rotation and control, and by which we can monitor any point 
\ on the pupil we wish is the Hartmann test. 

2. The Hartmann Test. 

a. General arrangement. 

In the classical application of the Hartmann test a mask of pre- 
cisely located holes is placed in the exit pupil of the optical system. 
Each hole lets through a particular ray converging toward focus. Two 
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photographic plates are exposed, one in front of focus and the other be- 
hind. By knowing the distance between the plates and by reading out on a 
comparator the positions of the rays on the plates, a slope can be 
determined for each ray in the pupil. 

This classical form of the Hartmann test is obviously not suited 
for remote testing. However, as illustrated in Fig. 20, only a simple 
modification is required. Plane A indicates the optical system exit 
pupil. We position lens B so that its rear focal point lies at the sys- 
tem focus. At location C we have the reimaged optical system pupil in 
collimated light space. In this relayed space we can sample a particular 
pupil location conveniently by means of a moveable aperture. Any depar- 
ture from a perfect system wavefront will show up in the sampled ray as a 
small angular deviation from collimated light. The slope error in the 
system exit pupil at this corresponding point can be found easily by 
multiplying the measured angular departure by the ratio of the relay lens 
focal length to the system focal length. This technique allows us to 
sample any pupil location we so desire using the full energy spectrum 
available, with an increased sensitivity equal to the focal length ratio. 

i. Automatic readout. 

We are now left to readout the angular deviations of these vari- 
ous rays. A suggested automatic technique involves the use of an image 
dissector tube of the star tracker variety. The arrangement 

is indicated in Fig. 20. The second lens D focuses the particular ray 
bundle we are sampling at its rear focal plane at some position off-axis 
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A 



pupil 

Fig. 20 . Hartmann sensor configuration. In summary the advantages are: (1 ) a priori 
geometric knowledge, (2) alignment insensitivity, (3) simple and efficient data 
reduction, and (4) high efficiency light usage. 


V. 
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which depends on the slope error and the sampled pupil position. The 
image dissector measures this off-axis location. This position divided 
by the lens D focal length is the magnified slope error at the pupil 
position. 

'fhe dissector tube consists of a photocathode that has been 
deposited within a transparent window. An electrostatic focusing system 
refocuses the electron image produced by the sampled ray bundle within or 
in the vicinity of a metallic plate containing a circular aperture in its 
center. Those electrons getting through the aperture produce a current. 
An orthogonal set of magnetic coils are phased 90° apart with equal 
deflection amplitude thus producing a circular scan of the electron image 
about the metallic circular aperture. Synchronous amplification of the 
resulting signal is carried out about the scanning frequency. If the ray 
image is perfectly centered in its scan about the circular disc, then 
there is no synchronous error signal produced. However,, if the ray 
bundle is displaced, an A.C. error signal is produced with some particu- 
lar phasing depending upon the ray position. From this we extract an 
X, Y-axis error signal which can then be applied as an offset signal to 
the deflection coils until the error signals are nulled. Knowing the 
calibration on the offset deflection we obtain the ray position. 

ii. Simulated automatic readout. 

This scheme has been simulated in the laboratory using a Risley 
prism scanner, a moveable pinhole and standard PMT. This is illustrated 
in Fig. 21. Lens B focuses a particular ray bundle on or near the 
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pinhole. This diffraction image is spun about the pinhole by a Risley 
prism whose orientation has been adjusted to give a circular spin 
diameter equal to the pinhole diameter. A mechanical chopper driven by 
the prism motor at 30 cps gives a reference signal phased on the x axis. 
Behind the pinhole is an E.M.R. 9858 photomultiplier tube. As the 
image is moved around the pinhole any offset would produce 
a 30 hertz variation in the P^^^ output current. As illustrated in Fig. 

22 this 30 hertz signal is preamplified and converted into a voltage 
signal. The voltage signal is then fed into a synchronous amplifier 
which first filters and amplifies the incoming signal and then splits 
the signal into two identical components 90° apart in phase. Each 
component is connected to the input of a phase-sensitive demodulator 
where the signal is mixed with the reference voltage signal to produce 
sum and difference currents. A filter at the output of the phase-sensi- 
tive demodulator rejects the high frequency (sum) frequencies and passes 
the zero or DC difference component corresponding to the component of the 
signal spectixim at the reference frequency. This zero- frequency current 
is the error signal output. We have two of these, one in-phase, indi- 
cating the x-axis displacement, and a quadrature signal, indicating the 
y-axis displacement. By moving the collection pinhole we can first zero 
out the x-error and then the y-error signals, read the micrometers and 
calculate our wavefront slope errors. Figure 22 is a photograph of the 
Hartmann test wavefront simulator being evaluated by the aberrated 
telescope. 



Figure 22. Hartmann Test Simulator and Aberrated Telescope Wavefront Simulator. 
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b. Theoretical discussion. 

We have begun to investigate several theoretical problems con- 
cerning the Hartmann test. These include an examination of the effect of 
sampling aperture size in the exit pupil on measurement of primary aber- 
rations and detection sensitivity. It also includes the basis for deter- 
mining the optimum size scanning pinhole for subsequent S/N ratio 
calculations. 

i. Effects of sampling aperture on third-order aberrations. 

The Hartmann test is used to measure the slope of a ray in the 
exit pupil of an optical system. By integration one may obtain the wave- 
front of the optical system under test. 

Ideally, we would like to make use of information from one point 
in the pupil, but since we are concerned about energy and diffraction we 
need a finite aperture for each ray. The question then arises: How does 

the aperture size affect our measurement of the ray slope, since we are 
integrating over a finite area of wavefront? In the following discussion 
we develop quantitatively a wavefront error measurement expression in 
terms of relative Hartmann apertinre size and aberration type. Diffrac- 
tion effects are considered in ii. 

The assumption is made that we haoe a measuring scheme that will 
locate the centroid of the image of the Hartmann aperture. Thus we are 
determining the average slope of rays in a Hartmann aperture and com- 
paring it to the slope of the ray at the center of the Hartmann aperture. 
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Fig. 23. Telescope system exit pupil. 


The exit pupil of the system is illustrated in Fig. 23. In a 

normalized system s 1. The center of the Hartmann ap'erture is 

located at Our analysis here restricts R s 1 - D, where D is 

the radius of the Hartmann aperture. Obviously Y = cos$ + y; 

X = R sin$ + X. 
c 

If S^XY x-direction slope of the wavefront at X, Y, then 

the average x-direction slope in the aperture centered around R^, $ is: 


+D 


S^^ydxdy 


xX.Y. 
1 1 


-D -(D^-x^)^ 
ttD2 
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Likewise the average slope in the y-direction is; 


yX.Y. 

•^11 


+D (D2-x2)^ 

L ! (1)2x2)^ 


SyXY^xdy 


ttD2 


Now let us consider several aberration types, and their measure- 
ment errors introduced in the Hartmann test. 


Defocus 


w = 


Y2) 

^xXY 

9W 

■ 3X “ 

2 %20 

^yXY 

3W 

3Y 

^ *^020 


In our small displaced aperture D, 

W = "^020 + x^ + y2 + 2 R^^x sin$ + y cos$J) 


xXY 


33T = '^020 f2x + 2R^sin$) 


xXY 




(2x + 2R^sin?>)dxdy 


xX.Y. 
1 1 


irD2 


ttD2 
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The first term integrates to zero and we are left with 


^xX.Y. ^''^020^i 

1 1 


or 


^xX.Y. ^xX.Y. 

11 11 


or the average x-slope in the aperture D is equal to the actual x-slope 
at the center coordinate of D for the defocus aberration. Likewise, 


Vx.Y. 

■'ll 


= S 


yX.Y. ■ 

■'ll 


Thus defocus in the absence of other aberrations causes no shift in the 
Hartmann image centroid. 


Tilt 

The wavefront aberration for tilt is W = W , , ,Y + W , , ,X where 

ylll xlll 

Wyjjl is a tilt of wavefront in the Y-direction. 

®xX.Y. ' Jx ‘ «xin(*-slope at X.Y ) 

1 1 

^x^Y " ^ " '^xlll 

II 

C = _ U 7 

x^Y 7 tD 2 xlll* 
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Thus, 

S . = S 

^XY ^X.Y. 

1 1 

We see that in the absence of other aberrations a tilt causes no shift in 
the Hartmann image centroid. 


Astigmatism 

The wavefront for astigmatism is 


W = ' " 222 ''^ 


S = 2W.,.Y. 

>'X.Y. 1 

1 1 


^ ~ ^ 222 ^^^ * * 2yR^cos$) 


s„ = = 2W„„„R cos<I> = 2W,^^Y. + 2yW~^, 

y^^ Y 9y 222 c 222 1 • 222 

i* 


2yW 


222 


u 


XY 


^^''^222^i * 2yW222)‘ix'iy 

ttD^ 


The second integral is zero. Thus, 


S = 2W._.Y. = S 

^XY 222 1 y 

1 1 


Again astigmatism does not cause a centroid shift. 
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Spherical and Coma 

Thus far with defocus, tilt and astigmatism we have seen that the 
Hartmann centroid is not shifted due to the sampling aperture size. We 
will see that for spherical and coma third-order aberrations, that this 
is not the case. That is, the finite aperture in the Hartmann test in- 
troduces a certain error of predictable form and magnitude depending on 
the relative Hartmann aperture size. 


Spherical 


W 


S 

X.Y 
1 





4W- .-X.R 2 
040 1 c 


4Wn.oY.R 2. 

040 1 c 


1 1 

In the Hartmann aperture: 

" ' "040("c" ^ x2 

[x sin$ + y 

/ / 

= * 


+ 


4"i’^040 



+ 2R^[x sin*? + y cos4'])2 

y2 + x2 + 2R 
' c 

cos$]) • (2x + 2R^sin$) 

/ / 

x^ + 2 X.x 2 + 2Y.xy] 

+ (x 2 .+ y 2 ) + 2X^x + 2Y^y]. 
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All odd functions integrate to zero as indicated by the /. We are left 
with: 


jj x^dxdy 

S = 8W_,„X. + 4X.W..„R 2 

^i i 


e2-n j-D 

j • p^pdpdS 
h JO 


+ 4X.W„._ 

1 040 itD2 


= S + 2D2x.W... + W-,.X.d2, 

Y x^ Y i 040 040 1 

1 i ^i 1 


By symmetry: 


= S ■ + 2D2y.W^.„ + W„,„Y.d2 

y^ Y Y ^ 040 040 i 

i i ^i i 


^ _ /m _ 

s„ dx. + S dY. 


■M^surad -x, , J ^ “‘i 

0 11 0 11 


w = w 

Measured 040 


* ^“o 4 oH • [v^v] 


^Measured ^040 V * 
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We see that when we measure spherical aberration with a Hartmann 
array we get an additional temi in our measurement, which takes the form 
of a focusing error with coefficient However, since we have 

obtained the functional form of the associated error we can now apply 
this knowledge to a correction term on our defocus measurement and thus 
compensate for the finite Hartmann aperture size. 
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Integrating we obtain: 


S = S 
"XY 


+ W_,d2. 


Xyv Y 131 

i i 


Since 




W. 


Measured 


M 


S + I S dY 
^XY ^ ^XY 


"'Measured "'iSl^^M^ "'iSI^m'^^' 


We thus see that in the case of coma, a finite size Hartmann ap- 
erture will introduce an error which takes the form of a tilt aberration. 
This information can be used as a correction on tilt measurement. 


Conclusions 

We see from the preceding discussion that if we are measuring the 
Hartmann centroid the only third-order aberrations that introduce errors 
are spherical and coma, and that the amounts of error are extremely small. 
We could predict the amounts of error for these aberration measurements 
by using the error formulas developed. We can also take advantage of our 
knowledge of the aberrational form of error to eliminate it. If we were 
monitoring an active system and were looking for changes, and we measured 
a change in spherical aberration we would also be able to predict the 
defocus type error that would be associated with it and thus not confuse 
it with an actual defocus change. The same argument holds for the coma 


and tilt error. 
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ii. Optimal signal recording geometry. 

In order to calculate the signal generated by the slightly off- 
center diffraction image of a sampled bundle of rays we must know the 
nature of the diffraction pattern formed on the image dissector's face, 
the impulse response of the tube and the dimensions of the collecting 
aperture within the tube. The convolution of these three quantities will 
give us the normalized amount of energy passing into the dissector for a 
given diffraction image collecting aperture displacement. From this con- 
volution curve we can determine the change in signal current for a change 
in position of the diffraction pattern on the tube face due to a slope 
error. Thus we will be able to make signal predictions based upon star 
source characteristics., f/no. of sampled pupil area, tube scanning aper- 
ture, sizes and tube impulse response characteristics. 

Referring to Fig. 21 let D be the sampling aperture size, f the 
lens focal length and X the average wavelength of the collimated light. 
The intensity distribution in the image dissector collecting plane for a 
unit amplitude disturbance is 


I(r) = 


2n4 


ir^D 
16X2f^ 


2J f— 1 

f— •) 
^Xf-’J 


The impulse response of the dissector tube is approximately Gaussian and 
can be written as 



G(r) 
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Here a is the Gaussian width of the response. The collecting aperture 
can be represented by H(r) = circ(r/£./2) , which is a function such that 

H(r) = 1 0 ^ T < j 

H(r) = 0 ^ ^ 1 

The energy getting through the aperture for a displaced diffraction image 
is 

L(r) = I(r)*G(r)*H(r), 

where * is convolution, and r is the coordinate of the diffraction 
image. 

Taking the Fourier transform (F) 

FL(r) = Fl(r) • FG(r) • FH(r) 

where 
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Multiplying and taking the inverse transform yields 


' ttD2£ 

L(r) = — j- 


J (2TrSr)J, (ir£S)e 
o 1 


-2-na^S^ 


cos 


‘ I - f (- 


dS 


where 


D 

^ A£ • 

We have evaluated this convolution for a series of possible c 
parameters, in particular for an image dissector with a Gaussian width 2o 
of 0.0035" and a collection aperture diameter £ of .020". This convolu- 
tion data was then used to determine the peak-to-peak change in signal 
per micron of displacement and finally the item of utmost importance: 
the change in signal per arc second deviation of the sampled ray, and a 
signal-to-noise figure of merit.* 

By plotting the S/N figure of merit versus focal- length (Fig. 24) 
we see that we do have an optimum choice of focal length. However, the 
function is slowly varying over a large focal length range. The real 
optimization choice should be made considering the displacement stability 
of the coils, which can be relaxed at larger focal lengths, and residual 
mechanical vibration problems which become worse at larger focal lengths. 
A good compromise for this tube is a focal length of about 50". 


*Signal current is defined as the peak-to-peak difference in the A.C. 
error signal generated by a ray displacement. Noise is determined by 
the total current generated by photons passing through the collection 
aperture Z. 
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TABLE III 


HARTMANN TEST SENSITIVITY VALUES 


eX 

D 

f 

% Signal Change 
per arc second 

S/N 

.01 

.1" 

10" 

4.3 

3.2 

.005 

.1" 

20" 

6.0 

4.1 

.002 

.1" 

50" 

7.2 

5.' 

.0015 

.1" 

66" 

8. 

5.2 

.001 

.1" 

100" 

10. 

4.1 

.0005 

.1" 

200" 

11. 

2.9 



Fig. 24. S/N figure of merit as a function of focal length (image dissector tube). 
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iii. Signal-to-noise for image dissector readout. 

Assigning we have optimized according to the above argument we can 
now calculate signal-to-noise ratios. 

For an extended red photocathode the total signal at the cathode, 
for a zero magnitude star assuming a 700 cm^ sampling area collecting 
aperture, is I = .50 x 10"^° amps on the photocathode, which corresponds 
to the maximum photocathode current the tube can handle. Now assuming the 
parameter optimization as described in the last section: 


(Shot Noise) 1^^ = /2eAfl 


Ij^ = /.5 X 10 -iO(2)(1.6 X 10-19)(1) 


= 4 X 10” amps. 


From above, for a 1 arc-sec detection, the signal is .061, 


Is = 30 X 10”^^ amp 


S/N = 30 X 10-^^ 

4 X lO-lS’ 


7.50 X 102 « 750, 


where dark current noise = 10“^^ amp for a .020” aperture. 

We compute signal-to-noise for other star magnitudes and summar- 
ize the results in Table IV. 
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TABLE IV 

BEST SIGNAL-TO-NOISE FOR DETECTING A 1 ARC-SEC SLOPE CHANGE 


Magnitude 

/ 

S/N 

Detection Certainty 

0 

750 

100% 

2.5 

250 

100% 

5 

60 

100% 

7.5 

6 

100% 


The above analysis has not considered noise problems associated with 
deflection coils and field effects in image tubes. 




CHAPTER IV 


FINE RANGE SENSORS 

1. The Zernike Test. - 

a. General discussion. 

The usual interpretation of the Zemike test is based upon the spa- 
tial filtering concept. Consider a disturbance in the pupil of 

an optical system as in Fig, 25. If the aberrations are very small then 
we can represent = 1 + i<|). The important observation is that 

the aberration information and pupil information are 90° out of phase. 

If, in addition, the aberrations are of a high spatial frequency type, 
the Fourier transform of the pupil function consists of a zero order, due 
to the pupil aperture, and a separated higher order spectrum 90° out of 
phase and separated from the zero order. (See Fig. 26.) We then argue 
that a phase delay of 90° introduced by a phase filter on the zero order 
will bring the higher order information in phase with the zero order. 

Lens C then re-transforms the filtered spectrum with the aberration 
information and pupil effects in phase so that the intensity pattern in 
the pupil contains a modulation proportional to the phase error in the 
pupil for small aberrations. Furthermore, by attenuating the zero order, 
the modulation effects can be improved. 
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♦ • 6 a e » 




Zero order 
at a phase of 0° 


Higher order spatial 
frequencies at a phase of 90° 


Fig, 26 . Spatial frequency distribution. 
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In the situation where we are looking for low order aberra- 
tions, such as may occur in figuring errors, thermal deformation errors, 
or misalignments, the problem is somewhat different. In this case the 
magnitudes of the aberrations are not so small that only a first-order 
expansion is permissible and, secondly, the aberration information is no 
longer separated out from the zero order diffraction pattern. These 
considerations lead to a slightly different interpretation of the problem 
as well as a more general formulation of the design of the Zernike disc 
itself, in terms of its size, the optimum phase shift, and the optimum 
attenuation. Zernike did some work on the low order problem in terms 
of a Zernike polynomial representation of the aberrations. However, he 
considered only the first-order term in the ■ expansion and did not gener- 
alize the problem considering obscuration effects. He also assumed the 
ir/4 phase shift and did not consider optimization of the disc parameters. 
What follows is a general development of the effects of a Zernike disc on 
the reimaged pupil. We show the relationship between disc size, phase 
shifts, and attenuations, and how these parameters affect signal-to-noise 
ratio, fringe visibility and the intensity in the reimaged pupil with and 
without various types of aberrations. We will also discuss the fabrica- 
tion of the disc, effects of fabrication errors and tolerancing. 

b. General development of intensity equation for reimaged pupil. 

. In Fig. 27 plane A represents the exit pupil, with fractional 


obscuration e of an optical system with radius a and focal length f . The 
exit pupil contains a general wavefront represented by The 
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normalized coordinate p ranges from e to 1. The pupil amplitude is then 
represented by: 




where 


H(^) = 1 


K- < 1 

X 


= 0 


> 1 . 

X 


At plane B, the Fraunhofer plane, the amplitude of the diffrac- 
tion pattern P is: 

P = *^iPo> 

where 


1 2tt 


= K. 




0 0 


K, 


la^ 

Xf 


2ir a 

^ ~ X f ^ 


and r, \p are coordinates in plane B, 

In plane B we place a circular mask of radius b, which produces a 
phase shift a within b and contains an amplitude attenuation (or relative 
amplification factor) T outside b. Thus, after passing through the mask 
P, the amplitude diffraction pattern becomes: 
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p. = , KjTP^[i - n(i)] * 


By rearranging we obtain: 


= KjTPq 


Ce'"-T^Wf^)* 


This amplitude field then passes through lens B and at plane C (the con- 
jugate plane to A) the amplitude of the disturbance Pq^ is: 



K2P' 


where 


2tt 


P» = 


P'e 


-iwr cos(i|)-6’) 


rdrdij; 


0 0 


2TTa*p 

X f 



q 


2ira 

Xf 


We are assuming lens B has an aperture sufficiently large to accept all 
scattered rays. Integrating we obtain: 



4t72KiK„ 



Ce^°-T)KjKj 



where 



3 2w 

I p^^-lwr 
0 0 
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is in general a complex function X(p,^) 


= = Q(p)e 


ii|) 


Q(p,ij/) = /Re X2 + Im X2 


'J'(pj'l') = Arctan 


Im X(p,if)) 
Re X(p,tji) 


Thus ; 


Pq" = Pq* " ltf2 • X(P) 


We now let (e^*^-T) = Ae^^, thus. 


i2 ^ 


1 + - 2T COSO 


tany 


sing 

cosa-T 


Substituting we obtain 


= T P, 


^A Q(p)e 


i(ip+y) 


What we have here is an expression for a two-beam interferometer 
where TP^* is the test beam containing the pupil aberration information 
and the second terra represents a reference wave whose araplitude ^ x2f2 '~ 
and phase y is determined by the Zernike disc properties. The factor ip 
as will be seen represents the effects of the aberrations on the refer- 
ence wave. As it turns out, \p is practically constant so that its effect 
is essentially one of shifting the intensity- gram by a constant. We will 
come back to the evaluation of X(p) in the next section. 
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Multiplying by its complex conjugate yields the intensity of 
the reimaged pupil for a unit amplitude pupil function. 


I.(p,e) 



j2 


+ 


1Q2a2 2TQA 
16X‘+F#‘t " 4X2F#2 


sin(5+\i;-<i)) 


(holds for E < p < 1) 
e > p > 1 

Our development here gives us an added bonus in the intensity 
prediction outside the pupil. We will see later how we can possibly use 
thfs information in order to •determine '■the zero aberration profile for an 
imperfectly manufactured Zernike disc. The question now is: how to 

choose the values of T and a to optimize the Zernike test performance? 

If we desire to operate on the midpoint of a. fringe so 
we have maximum range (-± quarter-wave) we want to choose ^ = 0. If this 
is the case we then have y = l + -^/l = tt/2. Then from tany = ^ v » we 

haye the constraining condition that T = cosa, and thus A = 1-T^ = sin^a. 
Now we have to consider fringe visibility. The visibility 

_ / 2QAT \ . / 2 . Q^A2 \ ^ / 2Qsin2g \ 

I 4X2F#2 j * r leX'tF#'* j \ 4X2F#2 / 

' * 16XWT ^ ~ 16X^F#‘+ 


where 


C = Y - tt/2 


I(p,e) = 


A^Q^Cp,6) 

16X4F#4 
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By differentiating and setting equal to zero, we obtain the condition for 
maximum visibility: 

cos2o = - 

Thus, by picking a to satisfy this equation and the transmission to 
satisfy the equation = cos^a, we guarantee a fringe visibility of 
unity and an operating point that gives us maximum range. We see that 
the choice of a depends upon ’ which depends upon wavelength, 

system F number, and the size of the Zernike disc. As we will demon- 
strate later, the effect of a large Zernike disc is to yield a reimaged 
pupil that has large intensity graduations for a zero -aberration condi- 
tion. A large Zernike disc also tends to mix aberration information into 
the reference wave. Smaller discs tend to yield a uniform zero- aberra- 
tion pupil profile with little aberration mixing in the reference wave. 
However the transmission factor is small and more energy is wasted. In 
Table V we have calculated the optimum phase shift and transmission 
required for several Zernike disc radii, corresponding to about 1/10 to 
1/2 the size of the diffraction pattern of an f/12 system operating at a 
wavelength of 0.6 ym, making use of expressions for Q derived in the 
following sections. Considering signal- to-noise requirements, pupil pro- 
file, and aberration mixing, the disc corresponding to about 1/3 the size 
of the central diffraction lobe appears to be a good compromise. 




67 


TABLE V 

ZERNIKE DISC DESIGN PARAMETERS AT 0.6 ym 


Disc Size 

Phase Delay a 

T2 

2y 

• 228X 

2% 

3y 

.206X 

7% 

4y 

.189X 

13.8% 

5y 

.178X 

18.8% 

c. Derivation of the 

effects of the Zernike 

disc on the reimaged 


pupi‘1 ■ 


The effect of the Zernike disc on the reimaged pupil manifests 
itself through the term X(p,9) = Q(p, 9) . One can evaluate the 
integral 


X(p,e) = 


b 2ir 

H 

0 0 


-iwr cos(iJ;-6) j j, 
U(r,\|»)e '^rdrdip. 


where U(r, is the amplitude diffraction pattern, numerically using a 
good deal of computer time. We have however, taken the approach of 
developing analytical expressions for both U(r,t;)) and X(p,9). The work 
has provided accurate analytical expressions for U(r,iJj) for aberrations 
of up to .22X and accurate analytical expressions for X(p,0) for Zernike 
disc sizes up to 2/5 of an Airy disc size. The range of accuracy could 
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be improved by expanding the expressions developed to include higher 
orders. The derivations have included the effects introduced by a system 
obscuration. 

First we develop analytical expressions for the dif- 
fraction pattern amplitude U(r,i|)), which includes aberration effects and 
the effect of a circular obscuration. The derivations that follow are 
done in terms of special functions Ajq(v) called Lambda functions. The 
Lambda function is defined such that: 




2^N! 




where is an Nth order Bessel function. Using the recursion formula 

for Bessel's function one can show that: 


A^Cv) = 


N-1 


4N(N+1) ^N+1^^^ 


* /■P^ 9 -.N M r • \ 

A(^)(l-p^) p COS 9 (-1) 


M JI_ 
N+1 


(2n) 


M N+1 


(V) 


(A) 


where +-+ represents a two-dimensional Fourier transform, and 

V = 27rr^ 

r^ = /X2 + Y2 



Y a 
o 


Y 


Xf 



We can also show that: 
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dA^(v) 

dv 



N+1 


d^Aj^Cv) 

dv 2 


2N+2 


dv 2 


9 Y 


- j V i[^y 


« r,i 

^N+2^ ^ ■ N+1 


] 


■U 


2ttcos.P|^- 2 V "N^J— 


'm 


C2tt)2 ,, ^N+2^^^ 

■^4 (l-cos 2 ti<) 


. (?j}l (i*cos 2«[|^A^^2M - Vl''')]' 


CB) 


The approach is now as follows. The pupil function e^*^ is 
expanded to the second term = 1 + i<j> - (j>^, where 


= Wjjpcos<|) + * ''^ 06 ^^ Wj^^P^coscj) 


+ W 22 P^cos^(|). 


This is a power series representation of the phase in the pupil where 
is a tilt effect, W^2 is a defocus effect, W^^ and are spherical 
aberration effects, W^^ is a coma effect, and W22 is an astigmatism 
effect. This expansion is substituted into: 
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2ir 

■ ^i+^-i[vpcos(e-«]p^pjg 
0 


where 


e'*^ = 1 + i<[) - 4)^/2. 

After substitution we arrive at a series of transforms of terms of the 
form p^^^^cos^O, where N and M take on various positive values. We can 
show that 


2N+M 

p cos 6 


■= 1 - N(l-p2) + (l-p2)2 


NCN-1) CN-2)(l-p 

•3; 


ill. ...] 


M 

• cos 0p . 


By taking each term in the power expansion and putting it in to the above 
form, applying the Fourier transform relationship (A) and finally making 
use of the derivative relationships (B) we arrive at the following 
expression for diffraction amplitude in plane B. 


U = + U,, + + U, , + U.,., + U,. ^ 

0 11 02 04 13 22 Interactions 


where U 


Interactions 


represent cross terms resulting from the (})2 term in 


the pupil function where: 


Uq = n[Aj^(v) - £2 A-(£v)] --diffraction pattern of the 

obstructed aperture. 


K 


2ir 

X 



Defocus: 


Uq 2 = n(iKWQ2)[Aj(v) - |a2(v) - enA^Cev) - | A2(ev)] 




2 


- *2''’’ " 3 ''3'''’ 


c'*/ajCeV) - Aj(cv) * i Aj(2V)j 


Third- order spherical: 


U 


04 


n(iKWQ^)[Aj(v) - A^Cv) + I AjCv) 

- e^^AjCev) - A2(ev) + j A^Cev)^ 
K%4^[Ai(v) - 2A2CV) + I AjCv) - A^(v) 


n 

2 


+ 5 AjM - e>» 


^AjCeI') - 2AjCev) 


+ 3 AjCgv) - A^(ev) 
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Tilt: 


Astigmatism: 


= n(iKW^^) (iv cosi|))| -j A^(v) - - e'+A^(ev) 


11 


*1_ 

_4 


4 2’ 


>] 


- K^W 2 
2 11 


2(^l(v) - 2 ^2^^^ 


+ - A^fv) |cos 2 ij; 


- e 


^ “ A2CEv)^cos2i|jj 


^22 " ^^^’^''^22^ 


|[ai(v) - i A2CV)] -H |[aj(v) - A 2 (v)]cos 2,J; 




Al(ev) - jA 2 (ev)j + yjA^Cev) - A 2 (ev)j cos 2 i|; 


— K^W 2 
16 ^ 22 


A^(v) - SA^Cv) + 3Aj^(v) + cos 241 • 


[ 2 A 2 (v) - eA^Cv) + 4A^(v)J 


) \ 


+ cos4ij^ 2A-(v) - 3A«(v) + A 


iCv)] 


- e®sAj(£v) - SA^Cev) + 3 Aj^(ev) + cos 2 i|' 



t 


2A^(ev) - 6A^(ev) + 4A (ev) 


Coma: 


* cos4<,|2A3(£v) - jAjCsv) * A,(,v) 


tJi 3 = - iAj 


A-Cv) 


- c-6 


4 13 


^^(4 ' X2 


AjCv) - I A^Cv) + A„(v) - 


* cos 2 <,j^AjCv) - 2 A 3 (v) * I Ajm] 


- ^8 


e llAj(ev) - - A (ev) + A (sv) 


- cos2^/Aj(cv) - 2A2(ev) . | A^Cev) 
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^INTERACTION 


■nK 2 Wo 2 Wo 4 


A, (v) 


I 


A,(v) - 4 


A4M 


- e®fA^(£v) - -| A2Cev) + A^Cev) - j A^Cev) 


- IIK2Wjj^WQ2(iv cosij;) 


4 ^ 2 ^^^ " 12 ^ 3 ^^^ 




- - l 2 





- W02W22) 


Ai(v) 


- A.(v) 


A3CV) 

" - 3 - 


+ COs 2 il; 


3 ^ 3 ^^^ 

AjCv) - I A2M .-V 


A^Cev) 


e^lAj(ev) - A 2 (ev)'+ ^ 


) 


'v‘ 


AjCev) 


+ cos2\1)IAj(ev) - -j A2(ev) + ^ 


iK^nv cosiKW^jWq^ + '^ 02'^13^ 


A 2 (v) AjCv) A^(v) 

—7— - — T— + 


24 


- e" 


lA^iev) A^Cev) A^Cev) 


4 


6 


24 



iK^Wo4Wi3nv cosij; 
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^2^^^ ^3(v) A^Cv) A^Cv) 


8 40 


- e 


A^Cev) A^Cev) 
4 4 8 


A^(ev) 

"~40 


iK^nWjjWjj 

16 


-vA^Cv) + Zvh^Cv^ cos\l> + ^vA^Cv) - vA^Cv)j cos3iJ; 


- E^^C-vA^Cev) + SvA^Cev)^ cos<|j + ^vA2(ev) - vA^Cev)^ cos3ij; 


iK2nWj2"^22 

_ 


(3*2 (v) - 2 Aj(v)j t i A^Cv) V cos* t Ajtv) - i Aj(v) 


i %<'')) 


V cos3i|) 


- £' 


^3A2 (ev) - 2A^(evJ + A^Cev)^ v cosijj. 


(a^Cev) - j A^(ev) + j A4 (ev)^v cos3t]; 


The above expression enables us to determine the diffraction pattern 
amplitude in the presence of third-order aberrations less than quarter- 
wave in magnitude, where the effect of a circular obscuration has been 


included. 
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We now want to determine the analytical expression that describes 
the effect of the Zemike disc on the reimaged pupil. We do this by sub- 
stituting the previous derived expression for U(r,iJ;) into the integral 
expression for X(p,0). We then make use of the fact that the A functions 
can be expanded in a power series where 


"N-1 


(ev) 


1 - 


e^v^ 

4N 


16N2(i + 2) 


We apply the three terms of the expansion to the unaberrated pupil U^, 
and the quadratic approximation to the aberration terms. We also make 
use of the recursion relationship for Bessel functions; 




2N J^(v) 


and the integral relationships: 


2tt 

cosM0e"^^*^^°^ d0 = (-i)^2fr Jj^(vp)cos MiJ;, 

0 

and 

1 

p^^^ j^cvpjdp = 

0 

After a good deal of tedious manipulation and evaluation we obtain: 
X(p,0) = X^(p,0) + X^^(p,0), 


I 
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where X^(p,0) is the effect of Zernike disc on the unaberrated pupil and 
XAb(p>0) is the effect rendered on the aberrated wavefront. 


X^Cp>0) =, 2ir^h^ 


_ bfai/_2_ , , , , 

wb 8 I w2b2 2 wb 


- 2ir^e^b^ 


— 


wb 


J3(wb)' 

wb 




J^(wb) 4J2(wb) SJ^Cwb) 
wb w2b2 ^ w3b3 


’‘AbarrationO-®) = ‘“'i 


02 


JjW ^ / 2 , , 

2»b ■ 12 ^w2b2 


J^Cwb) 

wb 


i2Tr2b2KWg2E‘* 


JjCwb) b2g2q2 (2J2(wb) J^Cwb) 

2wb 12 \ w2b2 wb 


b2772K2Wg2^ 


Jltwb) ^2^2 / 2 , _ , 

3wb “ 16 \w2b2 i^b 


+ b2TT2K2e^WQ2^ 


•’ 1 ^"'^^ bVji/ 2 , , ■'st"'’’ 

3wb 16 \ Vj2b2 2 wb 


+ iKW 


04 


>Jl(wb) ^2 2 ( 2 


3wb 




J^Cwb) 
wb 
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i2Tr2b2™ ,e6 

04 


Jj (wb) 
3wb 


b2^2q2 /2J2(wb) 
12 \ w2b2 


J^Cwb) 

wb 


b2Tr2K2w„ 2 

04 


J^(wb) 

5wb 


b^q^ 

24 



J2(wb) 


J^(wb) 

wb 


04 


Jl(wb) b 2 p 2 g 2 [23^{yih') J^(wb) 

5wb 24 V w2b2 . wb 


iq27T^KWj jCOS0 


J2(wb) _ ^2 k5 


J2(wb) 


4wb 


- „b 


iq27r2Kcos0Wj^ 


2 ■^2^' 

4wb 48 Iw^b^ 3^ wb 


b2 I_ k2w 2 
2 ^ 11 


_ bV / 2 , , . , 

2wb 12 \ w^b^ 2^” ^ wb 


J _ (wb) \ 1,2 2 


Pj(wb) 

I 2wb 


b2g.2q2 / 2 

12 U2b2 : 


(wb) 


J^(wb) 
wb ) 


+ cos26 


b^e2q2 

24wb 


J2(wb) 
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+ iq2tT^KlVj^2COS0 


.3 ^ sV (_ 2 _ ^ ... 

6wb 64 \w2b2 wb 


iq2iT2KWj2COs9e^ 


'^2^"'^^ ,, , q2e2bS / 2 , 

6wb ^ 64 \w2b2 “ wb 


J2(wb) 


K2ir2Wj2^b2 


JlCwb] ^ 
4wb ■ 20 


2 J3(wb) 

w2b2 wb 


+ cos26 


b2q2j^(wb) 
40wb 


Ml(wb) b2q2^2/ 2 . , J3(wb)^ oJb2e2q2 

— 20 — (^S2b7''2f"'>) - -Sb— j * 


- 


4wb 


+ iiT2KW 


22 


b^JiCwb) 2 / 2 _ 


2wb 


12 \w2b2 '^2*-'^^^ 


wb 




- e 


fb^JlCwb) J3(wb) 

2wb ■ 12 Iw2b2 '’ 2 ^”^' ■ wb 


^ ' . cos2e »‘'''' = ''’3t»« 


24wb 
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8 ""’'^“ 22 "' 


b2j^(wb) 




3wb 


b^q 

16 




J^Cwb)' 

wb 


+ cos20 


■’ 3 ^'*'') 


- cos4e 


£V 


J2(wb) J^Cwb) 
w2b2 6wb 




1 - J (wb) 

I o'- ^ 


- - n, (wb) 

w3b3 1 '■ 


( 


b2jj(wb) 

3wb 


b4q2g2 

16 




J2(wb) 

wb 


+ cos26 


b4q2g2 

24wb 


JjCwb) 


cos49 




J2(wb) 

w2b2 


JjCwb) 

6wb 





^•JlCwb)! 


plus cross terms. 

This equation has been programmed and the data run to obtain 
X(p,0) for cases that show the effect of the Zernike disc radius on X(p). 
As stated before, this information was used to calculate a and T and is 
tabulated in Table V. These data were then used to obtain the no-aberra- 
tion profiles for the reimaged pupil which appear in Fig. 28. These 
vividly demonstrate the tradeoff between disc size and pupil uniformity. 
For an f#12 system operating at an effective wavelength of 0.6 ym, a 3 to 
4 ym Zemike radius appears to be a reasonable compromise. 



Pupil intensity 
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Normalized pupil coordinate 

Fig. 28. Reimaged pupil through Zernike disk (no aberration). 






Pupil intensity 




Pupil intensity 
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Normalized pupil coordinate 


Fig. 32 . Reimaged pupil under the following conditions: X = 0.60 jum, 
F/no. = 1 2.000, b = 3.000, e = 0.4000, W 02 = 0.200, a = 1 .3025 
rad., 7 = 0.2605 



Normalized pupil coordinate 


Fig.- 33 . Reimaged pupil under the following conditions; X = 0.60 pm, 
F/no. = 1 2.000, b = 3.000, e = 0.4000, cjo 4 = 0.010, a = 1 .3025 
rad., T = 0.2605 







.0 


Normalized pupil coordinate 

Fig . '35 . Reimaged pupil under the following conditions: X = 0.60 jum, 
F/no. = 1 2.000, b = 3.000, e = 0.4000, W 04 = 0.100, a = 1 .3025 
rad., T = 0.2605 



Pupil intensity 


0.25 



Normalized pupil coordinate 


Fig. .'36 . Reimaged pupil under the following conditions: X = 0.60 /xm, 
F/no. = 1 2.000, b = 3.000, e = 0.4000, Wq 4 = 0.200, o = 1 .3025 
rad., T = 0.2605 
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The effects of aberrations upon the X(p,0) function were next 
analyzed. Using X(p,0) = Q(p, 0)e^’^^‘^’ we find a very interesting and 
important result. Using (defocus) as an example, the contribution to 
the Q(p,0) term, which governs the reimaged pupil shading, is insignifi- 
cant for Wq 2 and even less significant for higher order aberrations. The 
effect of the aberration is only on the <jj(p,0) term and this effect is 
for all practical purposes constant across the pupil. This tells us that 
the intensity changes we observe in the pupil will be solely a function 
of 4>, the pupil aberration, except for a constant average over the refer- 
ence wave passing through .the Zernike disc, which affects each point in 
the pupil by the same constant phase shift. Indeed if we look at higher 
order terms we find the effects to be even less. We have included 
Figs. 29-36 (graphs) of reimaged pupil intensity as influenced by varying 
amounts of defocus and spherical aberration from .OIX to .2X. A compari- 
son of these graphs with the zero- aberration case (Fig. 28) points out 
the sensitivity of the test. 


d. Calibration of real discs. 

Let us examine the graphs of the reimaged pupil. Figs. 29-36, 
concentrating our attention on the intensity profile outside the pupil. 

(r < e; r > 1) This area contains, practically speaking, little of the 
effects of the pupil aberrations for small aberrations and is in fact the 
functional representation of case of Zernike disc 

manufacture we will not have perfect square-well formation of the phase 
delay or exact dimensions. By measuring the diffracted intensity within 
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e and outside the pupil we will be able to fit a curve to determine 
AQCp, 0), for our real Zernike disc. This information will then be used 
in the pupil intensity equation to calculate the zero aberration profile 
of the real disc. Changes in intensity due to phase errors can then be 
related to this baseline. 

e. Tolerancing of discs. 

As mentioned before, the ideal transmission and phase delay can 
be calculated from the equations developed for the Zernike disc. Depar- 
tures from the design values in the manufacturing process will result in 
two main effects. The first effect is a departure from maximum visibil- 
ity, and the second effect will cause a shift in the zero- aberration 
operating point away from the zero condition. Figure 37 shows the effect 
on the visibility function V as a function of phase delay with three 
different transmission values. It can be seen that the visibility is 
indeed a slowly varying function and that for gross errors of phase delay 
and transmission its change is not significant. 

The significant effect of fabrication errors is illustrated by 
Fig. 38. Here we can see the shift in the operating point (the point 
midway between the maximum and minimum of a fringe) due to departures in 
phase and transmission from the design values. For example, if the phase 
delay upon manufacture was found to be 0.16 waves (instead of about 0.21 
waves which corresponds to the 0.00 operating point) the operating point 
has moved 0.05 waves from the symmetrical position. The range in one 
direction has thus been reduced from 0.25X to 0.20X. The prime purpose 
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Phase delay (waves) 

Fig.. 37 . Tolerance study for Zernike manufacture 
X = 0.6, e = 0.4, b = 3 



Operating point (waves) 



Phase delay (waves) 

Fig, 38. Tolerance study for Zernike manufacture 
X = 0,6, e = 0,4, b = 3 fim 
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then of tight tolerancing will be to maintain a maximum range of opera- 
tion (see Fig, 39). The other factor that influences our phase and 
transmission tolerancing is bandwidth. The shift in the operating point 
due to a finite bandwidth is ± AX/2X^ at the edge of the bandwidth where 
X^ is the nominal design wavelength. In Table VI we have tabulated the 
effects of error in the manufacturing process on the range of operation 
for the Zernike disc. For the erosion process, with careful control, and 
a 20 nm bandwidth the range of operation will be ± 1/8X in the worst case. 
For testing a well-corrected system the operating range can be reduced 
from ± 0.25X to a smaller range and the bandwidth can then be increased. 

For a deposition process, assuming ± 1/lOX deposition control we 
expect ± 1/5X in the worst case. 

f. Signal-to-noise calculations. 

One can show that in making a measurement the statistical signal- 
to-noise ratio due to photon noise will be: 

r. /V /AQTNatB 

' ~ /2 

V, A, Q, T, X, F# are properties of the Zernike disc and the 

incident wavefront as defined before. 

N is the number of photons per sec per cm^ per A° incident on 
telescope entrance pupil. 

a is the pupil area sampled. 

t is the sampling time. 


B is the bandwidth. 




Table VI. 


Erosion range ± 1/8 \ 


Error budget for Zernike disc manufacture. 

Assign! 1/120X transmission = ± 2% 

±1/10X phase delay = ± 1/5 X erosion depth 
±1/60X bandwidth = 20 nm 


Deposition range ± 1/5 X 


± 1/120 X transmission = ± 2% 
± 1/20 X phase plate = ± 1/10 X 


±1/60X bandwidth = 20 nm 


For a Zemike disc designed for C = 0, 


A2 = 1 - t2 


n /X, _ v/T(l-T2)^QNatB 

^ Hf# 

If we make two measurements, of different phases, the percent 
change in signal will be equal to VA(j> for small phase changes A<j). The 
signal-to-noise (S/N) must be greater than l/VAij). From this we can 
obtain an expression for the number of photons required to detect a phase 
change A(}). 


N = (NatB) 


8X^F#^ 

A((.‘+v3qT(1-T2)^ 


Figure 40 shows the minimum number of photons required to obtain 
A(j) = .OIX, A(}) = .02X for various size Zernike radii. 

Assuming a collecting area a = 100 cm^, a bandwidth of 20 nm and 
an A^ type star, we have determined the integration times required for 
various star magnitudes. Figure 41 shows the star magnitude vs. integra- 
tion time per pupil point for A^ = .OIX for 3 ym and 4 ym radii Zernike 
discs, and Fig. 42 shows the same for A<j> = .02X. These calculations have 
assumed a perfect detector. For a real detector with a quantum effi- 
ciency of about 0.1, the photon requirement would have to be increased by 
a factor of 10 and the integration times would also have to be increased 




Number of photons = , (NatB) 




CO 

E 

L. 

CO 

4 -* 

C/5 


0.1 0.33 


Time (sec) 

Fig- 41. Required star magnitude for X/100 detection 
A = 100 cm^, AX = 200 A 
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by a factor of 10. Thus for reasonable collection times in a real system 
we are talking about star magnitudes of 0 or less for the Zemike test to 
resolve Acj) = .OIX over the full ± 0.25X range. 

g. Fabrication process for Zernike discs. 

Below is a step-by-step procedure for manufacturing a Zemike 
disc, a procedure that we developed in our clean area solid state cir- 
cuitry facility. A good deal of the process is extremely sensitive to 
small dust particles, chemical contamination and environmental controls. 
These factors all influence the success or failure of the process. In 
summary dust must be minimized, care must be taken to avoid chemical con- 
tamination, and the temperature and humidity must be regulated in order 
to achieve repeatable results. 

i. Zemike Manufacture (Erosion Method) 

1. Clean substrate. 

2. Coat with aluminiim to specified transmission. 

3. Clean substrate. 

4. Spin on photoresist. 

Shipley: AZ-1350 

Spinner: - 5,000 rpra for 20 sec. 

5. Bake at 75°C for 10 min. 

6. Expose with mask in UV. 

7. Develop in AZ-Developer--45 sec. 

A. (1:1 dilution). 

B. Rinse in deionized water. 
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8. Bake 20 min. at 120°C. 

9 . Et ch . 

A. Aluminijm Etch. 

4 ml HNOj 

18 ml deionized H 2 O 
80 ml H^PO^ 


May heat to increase 
etch rate. 


B. Glass Etch. Rate (23*^0): 1.4 p/min. 


5:5:1 

H„0:NH,F.,:HF 
2 4 2 


10. Photoresist removal. 

J-100 at ~80°C 
Rinse in deionized H 2 O 
Dry - N 2 


11. Inspect. 
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ii. Zernike Manufacture - Metal Phase Shift Effect. 

We have developed the equations which describe the optimum transmission 

and phase delays. Ideally we would like the phase delay to be introduced 

by the erosion process and the attenuation introduced by the aluminum 

coating. However, because of the complex index of refraction and finite 

thickness of the metal we must make corrections to account for the phase 

delay introduced by the metal. We have calculated the phase advance 0, 

introduced by the metal for various values of thickness d and also the 

transmissions for various values of d. This information appears in Fig. 

2 

43. By referring to this graph and table V for T and o we can calculate 
the required depth of erosion in our substrate using the equation: 

t = (0 + Y ^ (n - 1) 

where 0 = phase advance produced by A1 in waves 

d = thickness of the A1 film 
\ = wave length 

0 = required retardation in waves 

n = index of refraction of substrate. 

As an example for a 3 y radius Zernike disc at 600 nm, we require that 
7 %, and a = .206 A. From Fig. 43, we obtain 0 = .121 X and 

11. nm. Taking n = 1.516 we obtain: 

t = C.121 + + .206) V (.516) 

t = .67 waves (at 600.00 nm) 
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used in Zernike manufacture. 
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2. Polarization Interferometers . 

a. General discussion . 

Polarization interferometers have generally been used up to 
now as devices for creating fringes. IVe can gain maximum sensitivity 
by taking phase information and converting it into states of linear 
polarization using a sensor which can read out the polarization state 
and thus obtain a phase reading. An instrument with a polarization 
readout in a Twyman-Green configuration has been built at the Optical 
Sciences Center and has demonstrated sensitivities of 1/500X. This 
is to be compared with a 1/50X sensitivity of the traditional two-beam 
fringe reading Twyman-Green interferometer. The Twyman-Green instruments, 
of course, operate with an artificial source which generates a nearly 
perfect reference wave. Since we are using a real star source we have 
to generate our own reference. Two possibilities are: 1) the use of 

some sort of polarization pupil shearing technique, such as the lateral 
shear, which we have chosen for laboratory purposes or, 2) the use of a 
method similar to the Zernike phase test, except operating in a polarizing 
sense, which can generate its own nearly perfect reference wave. The 
first concept we call the shearing polarization interferometer and the 
second the Zernike polarization disc interferometer. 

b. The Shearing Polarization Test . 

There are many types of shear (lateral, radial, inverting, etc.), 
which can be used to produce interference effects between two points in 
the pupil of an optical system. For laboratory simplicity we have chosen 
to illustrate, analyze and construct a polarizing lateral shearing device. 
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The basic idea is as follows: In Fig. 44 we have the exit pupil of 

the system which in gneral will be in some partially polarized state. 

We then reimaged this pupil at C through the polarizing interferometer 
as follows. By means of a polarizer we put the wavefront into a linear 
polarized state. This linear polarized pupil is then operated on by 
a series of optical elements to produce at C two images of the system 
exit pupil slightly displaced from one another such that one image is 
right-circularly polarized and the other image is left-circularly 
polarized. If we now look at a point in the sheared pupil image we 
have light right-circularly polarized from point A interfering with light 
left-circularly polarized from point B. If point A has phase and 
point B has phase Ag > the two oppositely-circularly polarized coherent 
beams combine to produce a linearly polarized beam such that 


where ijj is the direction of linearly polarized light. This is illus- 
trated in Fig. 45. We now use some sort of device to detect this linear 
polarization state and thus can obtain a direct measurement of the 
sheared phase (A^-Ag) across the pupil. 

i. Ways of Producing the Sheared Polarization States . 

There are many possible ways of obtaining the desired orthogonal 
circularly polarized sheared pupil condition at point C, and several of 
these illustrated below, with some brief comments. 
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Figure 46 shows a neat, compact unit where the shear in the 
pupil is produced by displacement of the two halves of the modified 
pentaprism along the polarizing beamsplitter plane A-B. Plane 
polarized light at 45° is then passed through the prism and sheared 
into two orthogonal states of polarization. These two orthogonal 
states then emerge at D where they pass through an achromatic quarter 
wave plate at 45° which in turn yields two oppositely-circularly 
polarized pupils at plane C. We have a compact unit where the shear is 
easily varied from zero to any value by movement along A-B. We require 
manufacture using very low birefringence glass so as to maintain linear 
polarization states. For greater stability one could replace the second 
right angle prism with another prism. The above arrangement is achromatic 
in collimated light. 

Figure 47 shows a compact way of obtaining circularly-polarized 
light by using two Wollaston prisms. The amount of shear is controlled 
by the separation of the two matched prisms with the scaling factor set 
by the prism angle 9, the unit being naturally achromatic with the main 
advantage in the compactness. Two disadvantages are in not being able to 
go to zero shear, and the large separation D required for large shears, 
which may result in uncompensated path differences and coherence length 
problems. This device appears most suitable for a small fixed shear mode 
where the two prisms are made into a single unit. In this case we have a 
single compact shearing mechanism. 
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In Fig. 48 we have a polarizing beamsplitter A-B. The shear is 
introduced by tilting a right angle prism about axis 0. The shear is 
given by 


S = 2t sin u [1 - ( 


1-sin^u 


> 2 _ 


sin-^u 


•) ] 


where t is the total thickness of the upper prism and u is the tilt 
angle. And 

S = 2tu [2^] 
for small values of u. 

In this case we can vary our shear from zero upwards. However, 
at large shears the variation is very nonlinear and wavelength dependent 
shear effects can be large. Low birefringence glass is required. 

In the simple arrangement shoivn in Fig. 49. a Wollaston prism is 
placed at the system focus. We have a fixed shear determined by the prism 
angle 9. Lens B will affect the polarization orthogonality of the shear 
and the arrangement is very sensitive to positioning of the prism. The shear 
itself is no longer achromatic. 

Summary 

There are other ways of introducing orthogonally-polarized sheared 
pupils either with birefringent materials or by using polarizing beam- 
splitters with tiltable glass plates or mirrors. Of the concepts pre- 
sented, based on maximum versatility and simplicity, the first would be 
preferred. 
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Fig. 46. Pentaprism polarization shearing device. 


Achromatic 



Fig. 47. Double Wollaston prism polarization 
shearing device. 


/ 





no 



Fig. 48. Polarizing beamsplitter tiltable prism polarization shearing 
device. 


Wollaston prism quarter-waveplate 



Fig. 49. Single Wollaston prism polarization shearing device. 
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Based on simplicity with limited versatility the second would be pre- 
ferred. For laboratory purposes and since we could use off-the-shelf 
components we have designed our shearing unit about this second concept. 


c. Ways of Reading Out Polarization States . 

We discuss here two possible ways of reading out polarization 


states. 



Behind this we position a Wollaston prism and a decollimating lens D, as 
in Fig. 50. If the original zero aberration position of the Wollaston 
prism was at 45°, we rotate this prism until the two oppositely-polarized 
intensities are equal. We can measure these intensities by use of an 
image dissector tube with deflection coils, in a synchronous amplifica- 
tion scheme (Fig. 50). We would thus rotate the Wollaston toward a null 
signal and encode the rotation, where a rotation of 1/60 radian (~1°) 
corresponds to a phase measurement ~ X/190. 

We have chosen to simulate this scheme in the laboratory because 
of its simplicity. However, instead of using an image dissector we 
simulate the dissector with two counter-rotating prisms as illustrated in 
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Fig. 50. Single channel polarization readout with image dissector. 


SHEARING POLARIZATION INTERFEROMETER 


Aberr. 

itar 

images 



Null 

meier 


Fig. 51. Laboratory simulation of shearing polarization interferometer. 
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Fig. 51. We thus sequentially sample one of two beams at 28 cps and 
synchronously amplify the difference signal. The Wollaston prism is 
rotated until this signal is minimized and this rotation is recorded. 

ii. The two- channel device . 

If we use the optical arrangement illustrated in Fig. 52 we get 
difference signals in each channel. In channel one, with the quarter- 
wave plate fast axis located at +45", we obtain a sinusoidally modulated 
signal where AI^^ = ~ sin(Aj^-A2) and in channel two with the 

quarter-wave plate fast axis oriented at -45° we obtain AI^ = K cosCA^-A^) 
The circuits are balanced for equal gain. Then the signals are divided 
so that 



tan(Aj-A2) , 


and then the arctanCAIj^/AI^) is taken. All of this is done electronically 
without the need for movement of any components. As with the other 
scheme it is insensitive to source fluctuations. This type of readout 
has been utilized in a Twyman-Green configuration instrument constructed 
at OSC and phase errors of 1/500X have been measured. 

d. Discussion of the Polarization Intensity Equations Including 
Alignment and Fabrication Errors . 

Below we drive expressions for the intensities of the sheared 
pupil signals as a function of (1) unbalanced signal amplitudes between 
two polarization states, (2) angular positioning of components and 





Synchronous differential amp 



pupils 


Fig. 52. Two-channel shearing polarization interferometer. 






(3) departure from perfect quarter-wave plate delay. 


Let us assume that the polarizer-shearing combination produces a 
shear in the y direction, as illustrated. Fig, 53, where the amplitude of 
one polarization is 

iA(x,y-|-) 


and the other is 




iA(x,y+j) 


and where S is the total shear. 


y 



Fig. 53. Sheared pupil geometry. 

These signals then pass through the quarter-wave plate oriented 

at some angle 9, with respect to the x axis. (Ideally 0 = ±45°, 

depending on whether we desire sin or cos modulation, as discussed in the 

previous section.) The plate advances the phase by ^ + <S for the fast axis 

2 


over the slow axis. 
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The field passing through the fast axis is: 


i(j + <S) iA(x,y - y) iA(x.y + 

Eg = e [E^e sin0 + E^^ e cos8]. 


The field passing through the slow axis direction is: 



We now orient a Wollaston prism at an angle 4' with respect to the x axis. 
Light through this prism then produces one signal polarized in the ip 
direction and another polarized at <|/ + 90®. j-. 



V 
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The output through the Wollaston is thus 


S = W = 1Eq12cos2(0-^) + |Eg^gJ2sin2(0-4,) 


+ 2Re EgEg^gQ cos(0-tj;)sin(e-ii;) 


^<^+90° "" ^!^+90^i|/+90* 1^0+90 Pcos2(0-,|;) 


+ E |2sin2(0-i|)) 
0 


- 2Re EgEg^gg cos(e-i|;)sin(e-iJ;) . 


We now substitute for E., and E. and assume 6 is small so that 

0 D+yu 


= i - 6. 


This yields 


= ^[E^2gi^20 + E^j2cos 20+ ^ E^^ sin20] cos 2 (6-\jj) 

+ [Ej 2 cos 20 +■ E,^ 2 gin 20 _ E^^sin 20 ] s in 2 ( 0 -ij;) 

- ^ E,/sin(20-2i())sinA 

+ <5, |^[Ej^ E/J cosA - (E„2_e^ 2) . ^ . ^ . q i ^ - j sin (20-2iJ;) 


cosA 



118 


I^^qqo = [E„ 2 sin ^0 + E^^cos^e - E^E^^ sin20]cos^ 

+ [E^^^cos^0 + E^^sin^0 + E^^ sin20] sin^ (0-iJ)) 

+ Ej^Ej^sin(20-2iJ^)sinA 

- 6[E^E,,cosA - (E^^2_E^2)Sn^ ^^.^^26-2,1;) 

We now substract the two signals: 


%+90° ■ \ = 2E^E^^sinAsin(20-2iJ;) 

- 2E_^E^^sin28cosAcos[20-2i/^) 

- 6S[E^E^,cosA - CE^^2 _e^ 2) ] sin (20-2i|;) 


i. Possible 

Operating Conditions 

for Two 

Channel Device 

Neglecting 6 

for the present : 



If 0 = 

( 

45°, = 0 

0 = 

0, ij; = 45° 

AI = 

2^ E^^ sinA 

AI = 

2^ E^/sinA 

If 0 = 

-45°, > = -45° 

0 = 

-45°, i\> = -45° 

• AI = 

2^E^^cosA 

AI = 

2^ E^^cosA 


We see we have two possible operating conditions for the two-channel 


cosA 


device. 
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ii. Departure from perfect quarter wave plate. 

In the one-channel device <{) is rotated until AI = 0, or 

" jniQ ■ 20+2i|;) + 6, - 2 ^E„cosA ’ : 

With the orientation of the quarter wave plate fast axis at 
0 = 45°, as far as tolerancing on 6 goes, we see that for small 

phase errors the second error term will be a small fraction of the 
first. We can manufacture a quarter ivave plate with 6 = ± 1/100 or 
less over 1000 A° bandwidth. 

• — 1 

. . 6 << -r^rTT and tanA = tan2ii + 6. 

ave 100 ^ 

The effect on small phase errors A = X/200 will be to introduce 
a constant offset error 6 , the constant offset being less than ± X/500. 
This constant correction would hold' beyond A becoming somewhat greater 
1/20X. Beyond this, the correction would get smaller for larger aberra- 
tions. If it were desired one could measure 6 over the passband of the 
instrument and determine the average 6, 6, and compensate for measured 
values . 

iii. Alignment of quarter wave plate . 

Now let us examine small errors in 0, A9. 0 = 45° ± A0. 

tanA = /n n ' "o A qV tan(90 - 20+2i|^) 

cos(9O-2A0) ^ 

1 

cos (2A0) 


tanA 


tan (2ip-2A6 ) . 
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The cos(2A0) effect is negligible for small A0 of a few degrees, 

. . tanA = tan (2i|^-2A0) 

A = 2i|; - 2A0. 

We see that the effect of the quarter-wave plate misalignment is to 
introduce a constant bias of 2A0 in the measurements , if A0 is kept to 
± 2 °. 


e. Bandwidth Considerations . 
i. General discussion . 

If 0 = 45° we arrive at an error signal AI = I (X)sih(A-2ii^) , 


A 

4 » —X 

x is physical path difference between 
the sheared points. 

4 = Kx 

K = 2ir/X 


AI = I (K)sin(Kx-2i(/)AK for monochromatic radiation of bandwidth AK. 


If we have a bandwidth of ± AK/2 and if we take into consideration 
the weighted spectral response of the detector, the shape of the bandpass 
filter and the distribution of the radiation we obtain an integrated 
effect. If the weighted function is symmetric about a we will show 

below that 


AI = sin(Kj^x - 2\|;)F(x, B(K),S, R)AK 


•Sv 


is the weighted average wave number. 


F(x, AK, S, R) is a function depending on the spectral 
distribution S, responsivity R, bandpass 
filter function B(K), and x. 
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It is desirable to have F =» 1 and a much slower varying function 
of X than the sine function. If we consider the case of narrow band- 
widths S, R can be assumed constant. If we take AK to be a square bandpass 
we find that : 

AI = sinfK X - 2ij;)sinc xAKAK. 

^ ave 

If AX = ± 250°A, and for the range of x, AKx ; 0, thus: 

sine xAK I const. ~ I. 


And 


AI r sin (K X - 2i{»)AK, 
^ ave 


for a narrow band where 


K + K . 
max min 


ave 


ii. General Theory 



In the above sketch F(K) represents the spectral energy dis- 
tribution multiplied by the filter response multiplied by the detector 
response. The error signal for a setting ip of our readout prism is 
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AI 



FCK)sin(Kx - 2<Ji)dK. 


In the monochromatic case we have seen that AI = 0, if 2i)j - Kx 

We would like to see the effect of this general radiation function F(K) 
on our error signal AI. 


sin(Kx - 2x^ 


1 

2i 



-i2ij; 
e ^ 


-iKx 

e 



AI 


AI 


J 


F(K)e"'^^dK - 


F(K)e‘^’''^dK 


Imaginary part of e 


-i24; 


F(K)e^’^dK 


i(Kx-2,j,) ^ 2-^0 

AI = Im e f F(K)e dK 

-(Ki-Kq) 

Here we have used the fact that F(K) is real, and the shift theo 
rem of Fourier transforms, where is the weighted average 

Now F(K) can be written as the sum of an even and odd function. 
F(K) = EF(K) + 0F(K) 

iCK X-2i|;) r r . 

.'. AI = Ime ° (EFCK)e dK + 0F(K)e dK 

The Fourier transform of an even function is real and even, and the 
Fourier transform of an odd function is imaginary and odd. Thus we 
obtain: 


AI = sin(K^x-2i|))F(EF) + i cos (K^x-2i|^) F(0F) , 
where F denotes Fourier transform. 
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We note that for any even distribution function our error signal 
is of the same functional form as in the monochromatic case, where is 
the weighted average wave number. In the case of a nonsymmetric distri- 
bution, the odd part of the distribution will cause in general an inability 
to obtain the null condition at K^x - 2ip = 0. Then condition for nulling 
now becomes : 

tan(K^x-2 1 {)) = - > (relationship A) 

In general we can thus make the bandwidth quite large if we take the 
spectral distribution function F(K) and correct for the null condition. 

In particular we have our highest sensitivity requirement in 
the vicinity of x ; 0. In this region F(0F) I 0, and we have 2\p = 
the monochromatic condition, as a solution above, and F(EF) only lowers the S/N 

iii. Summary. 

In summary our polarization interferometer can be used as a 
"white light" interferometer under the following conditions: 

(1) In the near zero aberration area with no corrections on the 
monochromatic case provided we use K = K^. 

(2) In the larger aberration area, with no corrections, if the 
asymmetry of the distribution is small. For example, if 

iF(OF) < 

F(EF) “ 60 ’ 

we would have 1/200X error. We also assume F(FE) 1. 

(3) In the larger aberration area if the spectral characteristics 
are known and we make corrections as provided by relationship 
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A, or if we are willing to accept lower accuracy for the 
larger aberration measurements, where the above relationship 
determines the measurement degradation if we know our dis- 
tribution and sensitivity functions. 
fi Signal-to-Noise Ratio 

In order to resolve a small change in signal we want the S/N 
ratio to be greater than the ratio of signal to change in signal, pro- 
duced by a phase difference, I/AI. For a very high degree of confidence 
99%, three sigma,, we require that S/N > 3I/AI. 

For A small, ijj Z 0, 6 = 45° 


AI : lA 



The 

approximate 


signal-to-noise results from taking tu)o measurements on an 
equal number of photons. 


S 

N 


»^#Photons 

/2 


DQE 


#Photons 

= NatB-E 


N 

= number incident per unit time per A° per 

sec. ~ 1000 

a 

= collecting area ; 100 cm^ 


B 

= bandwidth ; 500 A° 


E 

= polarizer effect = .5 (liriear polarizer 

X .5 (Wollaston prism) = .25. 

efficiency) 

DQE 

= Detective Quantum Efficiency of Detector 

s .33 



Required number of incident photons 
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Thus 


#Photons Required 



This relationship is graphed in Fig. 54. 


99% confidence 



Sensitivity 


Fig. 54. Relationship between the required number 
of incident photons versus sensitivity. 
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Using ; 

a = 

100 cm^ 


B = 

500 A° 


N = 

1060 photons/cm^sec A° for A^ zero magnitude star. 


We obtain a graph of sampling time versus star magnitude for various 
sensitivity requirements in Fig. 55. 



Fig. 55. Sampling times and stellar magnitudes required for three sensitivities. 


127 


g. The Effect of Shearing on the Low Order Aberrations. 

In the shearing interferometer we measure the difference in optical 
path between two points in the pupil. Below we derive expressions in 
terms of the shear S and the aberration coefficients. From these expres- 
sions we will be able to predict the minimum amounts of low order aber- 
rations we are sensitive to as a function of shear S, and instrument 
sensitivity. 


• > 




V 
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In the nonsheared pupil, if coma is oriented at angle 0, 
astigmatism at angle and tilt at angle ct, all measured from the y 
axis, the general wavefront expression is: 

W(x,y) = W 2 Q (x2+y2) + + 2x2y2 + y**) 

+ Wjj^x sina + Wj^^^y cosa + W^j^Cx^+y^Cx sin0 + y cos0) 

+ W 22 [^(x^+y^) + J 2 (y^-x^)cos 2 iJj + xy sin2\j)] . 

This is a normalized expression where x^+y^ <1. We now shear the pupils 
by an amount S, where S is a number from 0 to 2 (where S = 2 is a shear 
of a pupil diameter). This is illustrated above. At some point A, 
we now have two waves, W(+) and W(-) interfering. 

WC±) = W(x ± I , y) 

By substituting in the previous expression for W we obtain: 
w(±) = W 2 p{Cx 2 +y 2 ) ± -Sx + 

+ W^q{(x^ + 2x^y^ + y'^) ± 2x(x^+y^)S + j(3x^+y^)S^ 

±ixs3.^sn 

+ Wjj{x sin0 ± y s sina + y cosa) 
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+ W 2 j^{(x^+y^) (x sin8 + y cos0) ± y[(3x^+y^) (sin9) 

+ 2xy cos0]S 

+ -^[Sx sin0 + y cos0]S^ ± -^S^sin0} 

+ W 22 { [^^x2+y2) + 1 (y2-x2) (cos2ilj) + (xy sin2ii))] 

± y[ (l-cos2i|))x + (sin2i|))y]S 

+ |-[(1-cos2i|))]s2}. 

At point A, we measure the difference between W(+) and W(-). 

W(+) - W(-) = 2Wq2qXS + 4WQ4QxCx2+y2)S + W^^qXSS 

sinct + [ (3x^+y^) sin0 + 2xy cos0]S 

+ ^ W2jSin0S^ ■+ W 22 [ (l-cos2i|;)x + sin2i|)y]S. 

i. Sensitivity Dilution. 

We can make use of the expression above to determine our ability 
to measure various amounts of the different aberration types as a function 
of 1) our ability to measure W(+) - W(-) and 2) the amount of shear S. 
This information is summarized in Table 7 and Table 8. We see that in 
general the result of shearing is to reduce the sensitivity in measuring 
small amounts of aberration. This dilution, of course, is dependent upon 
the functional form of the aberration, and the amount of shear. For the 
LST a shear of . 2 -> .4 (the value of S) would be appropriate since 
it would give information over most all of the pupil. The important 
thing to recognize here is that the very fact of shearing dilutes the 



130 


sensitivity of any readout one chooses and that in order to make up for 
this dilution one must look for more sensitive forms of readout, as in 
the case of the polarization type readout, which we are investigating. 


ii. Small Shear Limit. 

As we make S small we can take the division 


W(-t-) - W(-) 
S 


= WjjSina + 2 ]/^ 20 ^ ^ ^W^^xCx^+y^) 
+ [ (3x^+y^)sin0 + 2xy cos6] 

+ W22[C1 -cos2iJ;)x + sin2ij;y]. 


This is merely the expression for 6W/6x, the slope of the wave- 
front (x component) at x,y, A corresponding shear in the y direction 
would then give the y component slope 6W/6y. 
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SHEAR 

.2 

1 

1.6 

ABERRATION 




DEFOCUS 

1/lOX 

1/25X 

1/16X 

TILT 

1/5X 

1/25X 

1/40X 

SPHERICAL 

1/20X 

1/50X 

1/6X 

COMA 

1/20X 

1/50X 

1/6X 

ASTIGMATISM 

1/5X 

1/12X 

1/8X 


TABLE VII 



Lateral shear. 

1/50X detection capability. 

static, fringe position reading. 


SHEAR 

.2 

1 

1.6 

ERROR 




DEFOCUS 

1/60X 

1/150X 

1/90X 

TILT 

1/30X 

1/150X 

1/300X 

SPHERICAL 

1/120X 

1/300X 

1/35X 

COMA 

1/120X 

1/300X 

1/35X 

ASTIGMATISM 

1/30X 

1/75X 

1/50X 


TABLE VIII 


Shearing polarization interferometer, assuming 1/300X capability with lateral 
shear (polarization or other phase reading scheme) . 
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3. Zernike Polarization Sensor (ZPS) 
a. General Discussion 

The ZPS represents a conceptual attempt to marry the Zernike test 
to a polarization readout. The result can be a sensor with at least the 
sensitivity of the Zernike test and a readout which provides for a 
direct phase measurement, in terms of linear polarization state. 

The experimental arrangement is indicated in Fig. 56. The wavefront 
from the system exit pupil first passes through a linear polarizer. At the 
system focus we place a birefringent crystal polished with its optic axis 
in the plane of the surface and adjusted to a multiple order thickness of a 
quarter wave. A small circular square well C~l/3 the Airy disc diameter) 
is eroded to a depth corresponding to a half-wave retardation between the 
ordinary and extraordinary waves. The area surrounding the square well is 
partially aluminized so that, with the Airy pattern centered across the well, 
the energy transmitted through the well is equal to the energy transmitted 
through the surround. As a result of this arrangment, we generate a refer- 
ence wave of one circular polarization sense from the inner part of the 
central core of the diffraction pattern, and a test wave of the opposite 
circular polarization sense from the outer part of the diffraction pattern. 
The lens following re-images the system pupil at plane C. In this plane, 
the two oppositely circular polarized waves interfere. According to the 
equations developed for the Zernike test, which are again applicable here, 
the test wave contains the wavefront deformation from each point on the 
pupil, and the reference wave is of essentially constant phase across the 
pupil, the constant dependent on the amount and type of aberration, and the 
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size of the eroded well. These two waves combine to produce a. linearly 
polarized state (the polarization direction depending on the phase 
difference at the point), at each point in the pupil. We measure this 
polarization using a device of the type discussed in Chapter IV and 
illustrated again in Fig. 56. We thus have a test with the simplicity 
of the Zernike test where we read out phase directly. 

b. Design Parameters for ZPS 

The optimum design parameters for the ZPS are obtained from the 

same equations derived for use with the Zernike disc. The effects of the 

diameter of the half-wave plate will be manifested in the phase variations 

in the generated reference wave. For relatively small constant effects, 

we want to keep this diameter about 1/3 the Airy core diameter. The 

transmission attenuation is introduced only because of its effect on the 

visibility function, since we are not concerned with the operating point 

2 

using the ZPS. The optimum T is obtained from Table V. 

c. Manufacture 

The success of the ZPS depends upon being able to manufacture a 
small half-wave plate. In order to eliminate tunneling effects we decided 
to look at materials with a high birefringent index. Two we have chosen 
to work with are calcite and rutile. 

We have prepared several samples of both materials cut with their 
optical axis lying in the planes of the polished surfaces. We are now 
trying to discover a procedure which will yield a smooth controllable sur- 
face erosion for these materials. We have tried many different acids on 
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calcite and have discovered one combination which is giving promising results. 
We have tried many combinations of acids under different conditions for ru- 
tile, but have not had positive results. 

We are now in the process of obtaining etch rate data on the calcite. 
After we have determined this we will continue on with the procedure in 
much the same way as in Zernike disc manufacture if we can find a substi- 
tute for the Shipley photo-resist (Shipley does not hold up under the acid 
combination used for calcite erosion) . 


Chapter V 


SUMMARY 

1 ■ Wavefront Error Simulator 

We have designed a F/12 wavefront error simulator that will enable 
us to test several types of figure sensors under varying . conditions of 
stellar magnitude, wavelength range, and aberration type and magnitude. The 
most critical part of this work was the design of an aberration generator. 
The design had to meet the specifications for an instrument that could 
simulate the aberrational effects of misalignment as well as of various 
asymmetries and orientations of small figure errors. We thus designed an 
instrument that could simulate errors up to several waves and at the same 
time have an incremental capability of around 1/lOOX. This instrument has 
been optically and mechanically designed, manufactured, assembled, and align 
ed. 

2. The Hartmann Test 

We have considered several methods of reading out the position of 
the ray being sampled and have decided that the best approach to position 
determination is the circular scan of a given ray around a small collecting 
aperture. We have worked on a program that takes the diffraction image size 
sampling aperture size and the impulse response of the detector, if appli- 
cable (i.e., star tracker), so that the readout parameters may be optimized. 
A S/N study has been pursued so that sensitivity vs. stellar magnitude can 
be predicted. Sensitivity predictions, optimum readout predictions and 
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other data reduction corrections and methods have to be verified in the 
laboratory. We are now proceeding with the experimental evaluation of the 
Hartmann test. 

3. The Shearing Polarization Interferometer 

We have examined the general capability of the lateral shearing 
interferometer and have reached the conclusion that a conventional lateral 
shearing device will not have adequate sensitivity because of the dilution 
of sensitivity inherent in a lateral shear. However, by incorporating a 
polarization type readout into the interferometer we can increase the sen- 
sitivity of the instrument significantly. Thus, in the laboratory, the 
primary objective in the construction of a lateral shearing interferometer 
is the experimental analysis of the polarization type readout using a star 
source. We have thus designed and manufactured a white light lateral 
shearing polarization interferometer in order to evaluate the polarization 
readout scheme. We have done some S/N studies and subsequent predictions 
on sensitivities and required stellar magnitudes. A preliminary error 
analysis on the system alignment and component fabrication has also yielded 
results. We plan to align the interferometer and evaluate it using the 
wavefront error simulator. 

4. Zernike Test 

We have worked on the general theory of the Zernike test and have 
developed expressions that can be used in the design of an optimum Zernike 
disc. We have also developed the theory that shows the effect of the finite 
size of the Zernike disc on the reimaged pupil due to the mixing of aberra- 
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tion information in the reference wave. This theory has been programmed 
and error analysis, S/N analysis, and some simulated experiments have been 
run. We have also worked on the technique for manufacturing the Zernike 
disc using an erosion process. Further work has to be done in the data 
reduction procedure for the Zernike test both for the 1/4X range and for 
larger aberrations. Work will be done on handling the effects of nonperfect 
square well phase shift and other disc irregularities. Work will also 
continue on the fabrication process. 

5. Zernike Polarization Test 

This test marries the Zernike test to a polarization readout, the 
possibilities being an extension of the sensitivity and range of the Zernike 
test. Much of the theory derived for the Zernike test is applicable here. 
The task that lies ahead now is an experimental effort on the manufacture 
of the device. One method we are investigating now involves erosion in 
calcite. If the manufacturing technique can be developed it should be 
followed by an extensive experimental evaluation of aberrations. 


